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Abstract: We provide exact solutions to the cosmological matter perturbation equation 
in a homogeneous FLRW universe with a vacuum energy that can be parametrized by a 
constant equation of state parameter w and a very accurate approximation for the Ansatz 
w(a) = wo + w a (l — a). We compute the growth index 7 = log /(a)/ log O m (a), and 
its redshift dependence, using the exact and approximate solutions in terms of Legendre 
polynomials and show that it can be parametrized as 7(a) = 70 + 7a(l — a) in most cases. 
■£3 We then compare four different types of dark energy (DE) models: toACDM, DGP, f(R) 

and a LTB-large-void model, which have very different behaviors at z^X. This allows us 
to study the possibility to differentiate between different DE alternatives using wide and 
£N| deep surveys like Euclid, which will measure both photometric and spectroscopic redshifts 

for several hundreds of millions of galaxies up to redshift z ~ 2. We do a Fisher matrix 
^vq analysis for the prospects of differentiating among the different DE models in terms of the 

growth index, taken as a given function of redshift or with a principal component analysis, 
with a value for each redshift bin for a Euclid-like survey. We use as observables the 
complete and marginalized power spectrum of galaxies P(k) and the Weak Lensing (WL) 
power spectrum. We find that, using P(k), one can reach (2%, 5%) errors in (wo,w a ), and 
(4%, 12%) errors in (70,7a), while using WL we get errors at least twice as large. These 
estimates allow us to differentiate easily between DGP, f(R) models and ACDM, while it 
would be more difficult to distinguish the latter from a variable equation of state parameter 
or LTB models using only the growth index. 
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1 Introduction 

With the turn of the millenium we have entered a new era in which cosmological observa- 
tions have improved to the level that we can start to define a Standard Model of Cosmology 
based on the CDM paradigm plus some sort of vacuum energy responsible for the observed 
dimming of distant supernovae. The nature of either Dark Matter (DM) or Dark Energy 
(DE) remains a mystery, in spite of the improved determinations of their contribution to 
the total energy density of the Universe. 

While the nature of Dark Matter seems less uncertain (most cosmologists are in favor 
of a particle physics origin), that of Dark Energy is still unexplored territory. In the last 
decade there has been a plethora of proposals to account for the observed acceleration of the 
universe. All these proposals fall into four main categories: i) the inclusion of some extra 
field (scalar, vector or tensor), coupled or not to the rest of matter, like in quintessence, 
chameleon, vector dark energy or massive gravity; ii) the extension of GR by inclusion of 
higher order terms in the Einsten-Hilbert action, like f(R) theories, Gauss-Bonnet terms, 
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etc.; iii) the modification of gravity on large scales by introduction of extra dimensions, like 
in the Dvali-Gabadadze-Porrati model, Kaluza-Klein gravity, etc.; iv) the reinterpretation 
in terms of a nontrivial spatial geometry, like in large-void inhomogeneous LTB models. 
For a recent review on DE modeling, see for example [1, 2]. 

All of these proposals have very specific predictions for the background evolution of 
cosmological space-times, and most of them can be well fitted to the present observations 
with just a few phenomenological parameters: the equation of state, the speed of sound, 
the coupling between DM and DE, bulk viscosity, etc. However, in order to discriminate 
between the different alternatives it has been realized that one has to go beyond the back- 
ground evolution and start to consider also the theory of linear cosmological perturbations 
and parametrize their evolution in terms of the growth function and growth index, as well 
as the shift parameter. 

At the moment, the main observables used to constrain the dark sector are the cepheids 
and supernovae magnitudes for the determination of the expansion rate as a function of 
redshift, see [3] for the most recent measurements; the power spectrum of matter and 
the baryon acoustic oscillation (BAO) scale for the determination of the matter content 
as a function of redshift, see [4, 5], together with the cosmic microwave background for 
the determination of global spatial curvature and the asymptotic values of cosmological 
parameters [6], and the weak lensing shear power spectrum [7], or the ISW-galaxy cross- 
correlation [8], for consistency of the whole scenario. Most of these measurements are 
rather preliminary and suggest a detection of DE at the 2 to 3 sigma level. However, a 
significant improvement is expected in the present decade thanks to Planck [9, 10], DES [11], 
BOSS [12] and, in the future, the Euclid survey [13]. 

In this paper, we study the prospects that a survey like Euclid would have in distin- 
guishing between the four main classes of DE models. In the process, we find exact and 
approximate solutions for the growth index in terms of simple functions, for ACDM mod- 
els with a constant and variable equation of state parameter. We then propose a simple 
parametrization of the growth index that fits well the recent history, except for extreme 
models like f(R). We note that while the background parameters H(z), Qm(z) and w(z) 
seem to be rather similar to those of ACDM, the growth index can differ significantly for 
most classes of models. 

The parametrization of the density contrast with the growth index was first introduced 
by Peebles in 1980 [14] . We knew that the rate of growth of structures should be a function 
of the matter density. Several parametrization attempts including a power law expansion 
or simply the square root of the matter density parameter did not quite fit the data. 
Finally, the 7 parametrization was the most widely accepted one; however, we will have 
to wait for a direct measurements of the growth index, which will only be possible with 
the next generation of experiments. It is therefore tantalizing to explore the possibilities 
of distinguishing between different DE models with a better determination of the growth 
index than what we have at present. That is the reason why we study a survey like Euclid, 
that will allow us to obtain information not only about the matter distribution (power 
spectrum P(k) of perturbations) but also about the weak lensing (WL) spectrum. 

We perform a Fisher matrix analysis of the sensitivity of a Euclid-like survey to the 
growth index using the marginalized and complete power spectrum P(k) and WL observ- 
ables, and find that one can improve significantly its determination if we know the specific 
form of the redshift dependence of the growth index. We find that using P(k) alone, one 
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can reach (2%, 5%) errors in (wo,w a ), and (4%, 12%) errors in (70,7a), while using WL we 
get errors at least twice as large. These estimates allow us to differentiate easily between 
DGP, f(R) models and ACDM, while it would be more difficult to distinguish the latter 
from a variable equation of state parameter or LTB models using only the growth index. 



2 The background equations 

Here we review the basic equations for the relevant background quantities. The evolution of 
the dark energy can be expressed by the present dark energy density Q,de and its equation 
of state parameter: 
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The Hubble parameter, H(a), is 
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where the subscript denotes the present epoch, and we are assuming global spatial flatness 
{K = 0). The angular diameter distance becomes 
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The Eqs. (2.4) and (2.5) imply that we limit ourselves to a flat universe which is filled with 
matter (denoted by the subscript m) and a general dark enery fluid. 
The total matter and dark energy densities are: 
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3 The perturbation equations 

In this section, we consider linear perturbations about a spatially-flat background model, 
defined by the line element: 



ds 2 = a 2 (r) [- (1 + 2*) dr 2 + (1 - 2$) dxidx* 



(3.1) 
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where r is the conformal time. From Eq. (3.1) it is clear that we are working in the 
Conformal Newtonian (or longitudinal) gauge, which influences the evolution of the per- 
turbations on big scales (especially those larger than the Hubble horizon). However, on 
smaller scales the choice of the gauge is less important, and we expect the perturbations 
to evolve independently from the gauge choice. 

The perturbation equations for a general fluid with an equation of state parameter 
w = p/p are [15, 16]: 
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where 5 = Sp/p is the density contrast, V = ikfT^jp is the scalar velocity perturbation 
(see [17]) and the primes denote derivatives with respect to the scale factor a; also 5p and 
a are the pressure perturbation and the anisotropic stress, respectively, and they depend 
on the characteristics of the fluid itself. 

In this paper we are interested on the evolution of the matter density field, which 
implies that w = bp = a = 0, at small scales (but sufficiently large for linear perturbation 
theory to still hold). Furthermore, we also assume that the dark energy fluid enters only 
at the background level. This assumption is usually violated as dark energy perturbations 
may influence the evolution of the matter density field especially when the dark energy 
sound speed is small [18-20], since it is then able to cluster at small scales. However, 
under the assumption of a high value of the dark energy sound speed (as is the case for the 
quintessence model, where c 2 = 1) we expect this assumption still to hold because dark 
energy is able to cluster only on very large scales. 

Under these assumptions, Eqs. (3.2) and (3.3) read: 

(3 - 4) 

The above equations are linked together through the gravitational potential $ 

fc 2 d> = -47rGa 2 p m [b m + ^V^\ (3.6) 

and \I> = <E>. 

Joining Eq. (3.4) and Eq. (3.5), the master equation for the matter denstiy contrast 
becomes: 

3 

a 2 (5^(a) + (3 - e(a))a5' m (a) - -VL m (a)5 m (a) = , (3.7) 

where e(a) = — d log H{a)/d log a. The exact growing mode solution of the above differential 
equation, for a constant dark energy equation of state parameter, w, is [21, 22] 
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where 2^1 are Hypergeometric functions, see Ref. [23], and the full solution, including de- 
caying modes, is reported in the Appendix A. Equation (3.8) can be further simplified 
using trasformation formulae for the Hypergeometric functions [23] into 
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(3.9) 



It is worth having a further look at the structure of the solution found. Eq. (3.9) is 
composed of two terms: the first one is the usual scale factor a, which is the solution we 
expect if dark energy had been neglected also at background level; furthermore, Eq. (3.9) 
reduces to the classical solution in the matter domination era; this can be clearly seen if we 
set Q m (a) —> 1, then the Hypergeometric function is 1. The second term instead, contains 
all the information about the dark energy fluid: via a direct dependence on Q m (a) and the 
Hypergeometric function, whose dependence is not straightforward. It is worth noticing 
that the slowed contribution to the dark matter comes from the term f2 m (a) - 3^ but its 
contribution is too big and it suppresses too much the growth of matter perturbations, 
the Hypergeometric function in Eq. (3.9), instead, is a function that slowly increases with 
the scale factor, from 1 for a <C 1 to just 1.16 for a ~ 1. This compensates for the extra 
suppression of the f2 m (a)~3^ term. 

In order to evaluate the growth rate, we start using Eq. (3.9), and write it as, see 
Appendix A, 
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This function is simpler to evaluate for a given model, and has the explicit dependence 
of the growth factor as a function of O m (a). We will use it in the next subsection on the 
gamma parametrization. 

3.1 Varying equation of state parameter 

We can extend our discussion also to a varying dark energy equation of state parameter 
w(a). In this case there are no exact analytical solutions for the matter density contrast 
as for the case of a constant w; the main problem is that there is no direct transformation 
between the scale factor a and the new variable u, see Appendix A. Here we assume the 
equation of state parameter to be: 

w(a) = w + w a (1 - a) , (3.11) 
for which the matter density parameter can be integrated 
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We can also integrate the master equation in this case, and although we have found no 
exact analytical solution, we found an approximate solution which is within 0.1% of the 
numerical solution for the whole range of (wo,w a ) values. In this case the growth rate 
reads 

f(a) = (a) r .„ , 4 ■ (3.13) 



p5 / '6w(a)+w a a/6w 2 (a) 
"l/6w(o) 


^m 1/2 (a) 


j-,5/6w(a)+w a a/6w 2 (a) 
-l/6io(o) 


^m 1/2 (o) 



This is a very compact expression which will be very useful in the next subsection. 
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3.2 The gamma parametrization 

A potentially more promising place to look for the impact of dark energy is the way in 
which the matter perturbations grow over time. Often the impact of the dark energy on 
the growth rate of the matter perturbations is parametrized in terms of the growth index 
7, defined through [14]: 
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Such a parametrization assumed 7 to be constant and, moreover, with a value that was 
very approximately 7 ~ 0.6 for general relativity plus a cosmological constant A, i.e. for a 
constant equation of state parameter w = p/p = — 1. Furthermore, detailed studies went 
further and computed the 7 parameter for a constant but arbitrary w, [24]: 
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which reduces to 7 = 6/11 ~ 0.55 for w = —1. The last equation does not take into 
account the dependence of 7 on the matter density parameter. There is an expression for 
the growth index 7 which does depend on O m (a) < 1, to first order, see Ref. [25], 
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which reduces to the well known result (3.15) in the limit a — > 0. However, we have found 
in Eq. (3.10) an exact solution for all a. Using Eq. (3.14) we have 
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In particular, Eq. (3.17) depends not only on w, but also on O mj o- For the fiducial values 
w = — 1 and £l m ,o = 0.25, one finds 7(0 = 1) = 0.556 instead of 7(0 — > 0) = 0.545. For 
present day purposes, with galaxy surveys providing at most a few percent accuracy on 
the growth parameter, this difference — of order 3% — may seem academic. However, for 
future surveys like PAU [26], LSST [27] or Euclid [13], where we will have tomographic 
reconstruction of the past history in both the matter distribution and the expansion rate, up 
to redshift z ~ 2, these differences may begin to play an important role as a discriminator 
between standard GR with a cosmological constant and e.g. modified gravity theories like 
f(R), or quintessence models. 

The same discussions apply when we want to evaluate the growth index when the 
equation of state parameter is a function of the scale factor, 
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In the next section we will discuss the growth index for different Dark Energy models. 
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4 The growth index for different Dark Energy models 



Here we discuss a range of models of dark energy which have very different behaviours 
for the growth index 7 as a function of redshift. The models we consider here are f(R) 
theories, the DGP model, LTB cosmologies and wCDM. With this choice of models we 
wish to cover the different possibilities to explain dark energy: modified gravity theories 
(f(R) and DGP models), modified geometry (LTB models) and modified matter through 
wCDM, which can be associated with, for example, a quintessence model through the 
correspondence between f((f>) and w(a). 

4.1 The Dvali-Gabadadze-Porrati model 

The brane-world scenario offers an alternative to explain the current acceleration of the 
universe. In the DGP model [28], gravity lives in a 4 dimensional brane with the rest 
of matter fields for scales A < r c while it seeps into a fifth dimensional bulk for A > r c , 
weakening its strength in our brane. For such models, the Hubble parameter is given by 
[29] 

H(a) = H y/n^+Jn rc +n mfi a~ 3 (4.1) 



where Q rc = l/(4r^i?o) = (1 — ^m,o) 2 /4. The Poisson equation is modified in these models 
as follows: 
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For the growth index, we have used the approximation found by 
7 + 5n M (a) + 70^(a) + 3Q^(a) 
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while the effective equation of state reads 
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4.2 The Starobinsky f(R) model 

The simplest way to modify general relativity is by substituting the Ricci scalar R with a 
function f(R), in the Einstein-Hilbert action, [32]: 
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The equivalent of the Friedmann equations for a perfect nonrelativistic matter fluid for this 
action are given by 

m 2 ' (4.7) 

Pm + F- HF, 



-2FH 



-7- 



where F = f'(R). The Ricci scalar R is expressed in terms of the Hubble parameter as 
R = 12H 2 + 6H for a flat background. It is convenient to introduce a set of dimensionless 
variables to study the dynamics of f(R) gravity [33]: 
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With these variables we obtain the following closed system of differential equations 
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The linear perturbation equation for the matter density contrast in the comoving gauge 
(where the matter velocity vanishes), in terms of the dimensionless variables, x\, X2 and 
x 3 , is [34] 
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where the differential equation for SF = 5F/F is 
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Solving the coupled differential equations (4.10), we can obtain the background functions 
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Then solving equations (4.12) and (4.13) numerically for a given scale k, we find a solution 
for 5 m , and compute from it the growth index 7. 

In this paper we study the Starobinsky model [35] where f(R) is given by 



f(R) = R + XRo 



1 + 



R^ 



(4.15) 



where A and n are two positive constants and Rq corresponds to the present value of the 
Ricci scalar. To be in agreement with observations, we take n = 2 and A = 2 [36]. We also 
take k = 0.16/i/Mpc in Eq. (4.12) to calculate the density contrast and the growth index. 
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4.3 The Lemaitre-Tolman-Bondi large-void model 

A plausible alternative to explain the current acceleration of the universe are inhomoge- 
neous universe models in which the acceleration we seem to perceive is caused by our 
position as observers inside an underdense region of space. One of the simplest models 
to study the effect of such large inhomogeneities is the spherically symmetric Lemaitre- 
Tolman-Bondi model [37-39]. In this model, the metric is given by 

ds 2 = -dt 2 + X 2 (r, t) dr 2 + A 2 (r, t) dti 2 , (4.16) 

where dQ 2 = d6 2 + sin 2 Odcf 2 and the equivalent of the scale factor now depends on the 
radial coordinate as well as time. We can find a relationship between X(r,t) and A(r, t) 
using the the — r component of the Einstein equations: X(r,t) = A'(r, t)/Wl — k(r) 
where a prime denotes a derivative with respect to coordinate r and k(r) is an arbitrary 
function that plays the role of the spatial curvature parameter. 

To find the growth index in LTB cosmologies we must study linear perturbation theory 
in inhomogeneous universes. Due to the loss of a degree of symmetry, the decomposition 
theorem does no longer hold. This means that, in general, our perturbations will no longer 
decouple into scalar, vector and tensor modes. A study of the perturbation equations in 
this scenario using a 1+1+2 decomposition of spacetime can be found in [40]. However, if 
the normalized shear e = (Ht — Hl) I (2-ffr + Hi) is small, as observations seem to confirm 
[41], we can use the ADM formalism and express our perturbed LTB metric as 

ds 2 = -(1 + 2$>)dt 2 + (1 - 2^) lij dx i dx j (4.17) 

where 7^ = diag{A 2 (r, t), A 2 (r, t), A 2 {r, t) sin 2 9}. Within this formalism, the evolution 
equation for a pressure-less fluid in the conformal Newtonian gauge (where the absence of 
anisotropic stresses gives 3> = \£) is given by 

$ + AH T $ + (4F T + 6ilf,)$ = (4.18) 

where we now have two different expansion rates Hx(r, t) = A/A and HL(r,t) = A' /A' 
which correspond respectively to the transverse and longitudinal expansion rates. The 
growing mode solution of equation (4.18) is given by 

*(r,t) = *o(r) 2 Fi[l,2 ) ^;u] ) (4.19) 

where u = k(r)A(r,t)/F(r) and F(r) = H 2 (r)QM{r)A 3 (r,to) specifies the local matter 
density today. With this solution we find the density contrast, 

5(r,t) = ^$-$(r,t) (4.20) 
r 

We can now calculate the growth index as before, noting that now the matter density 
parameter is a function of redshift via both time t and the radial coordinate r. In LTB 
models, this is in principle an arbitrary function which must be fixed in each case. In the 
case of the constrained GBH model [42, 43] the parameters are given by 

n M (r) = 1 + (n<? - 1} l-taiih[(r-ro)/2Ar] 

u ; 1 + tanh[r /Ar] v ' 
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with 



ro = 3.0 Gpc , Ar = 1.5 ro 
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(4.23) 



where these values have been chosen to best fit the supernovae data [41, 44]. Within this 
model, the growth rate, i.e. the logarithmic derivative of the density contrast, is given by 
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where Q m (z) is the fraction of matter density to critical density, as a function of redshift. 
This function (4.24) is identical to the instantaneous growth function of matter density in 
an open universe, where the local matter density £Im is given by Q m (z) at that redshift. 
This is a good approximation only in LTB models with small cosmic shear, see Ref. [45]. 
Alternatively, we can write the growth function in terms of Legendre polynomials, 



f(z) = 



p 



-5/2 
-1/2 



n- 1/2 ( z ) 



p, 



-5/2 



1/2 



n m l/2 (z) 



(4.25) 



which can also be written in terms of ordinary functions, 



/(*) 



9tt(l 



u 



+ 6(1 + 2u 2 )Vu 2 - 1 arcsinh v / (n - l)/2 



2u(u 4 + u 2 - 1) 



I2u 2 y/v? 



1 arcsinh-^/ (u — l)/2 



(4.26) 



n m 1/2 (z) > 1. 



where u 

The equation of state w(z) has been obtained for this model from the expression 

diogCn-Ho)-!)- 1 



w(a) 



d log a 3 



a£l' m (a)/n m (a) 
3(l-O m (a)) 



(4.27) 



which was used in Fig. 1. Note that the rate of expansion H(z) for this model is similar 
to that of ACDM, which explains why it fits the SNIa data [46, 47]. 

4.4 Comparison of the growth index for the different models 

Now we are ready to compare the results of the different models studied here. In Fig. 1 
we can see the Hubble rate, the matter density parameter, the equation of state and the 
growth index for all the models studied. In the case of wCDM we have taken wq = —0.9 
and w a = 0.2. As it is seen, even though the expansion history, matter density parameter 
and equation of state are quite different for the different models considered, most of them 
have a similar growth index, being the biggest difference between DGP, f(R) and the rest. 

This difference in growth indices comes from how linear perturbations in the matter 
density are described in each of the models. In the case of f(R), if we neglect the oscillation 
mode of 5F relative to the mode induced by matter perturbations 6 m , we can obtain the 
following approximate equation for matter perturbations: 



5 m + 2H8 m - 47rG e sp m 5 m ~ 



(4.28) 



1 The matter density in LTB model is given by p(r,t) = F'(r)/A'(r,t)A 2 (r,t). Note that this is different 
from Sl M (r) = F(r)/A 3 {r,t )Hi{r), which gives the mass radial function today, see Ref. [42]. 
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(a) Hubble rate 



(b) Matter density 
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Figure 1: Hubble rate, matter density parameter, equation of state and growth index as 
a function of redshift for different DE models. The solid lines represent the ACDM model, 
the short dashed lines the wCDM one, the dotted lines the DGP model, the dotted dashed 
lines the f(R) theory and the long dashed lines, the LTB model. 



where 



k 2 F' 



= ^ I — -Mr (4-29) 

where here a prime means a derivative with respect to R, and G is Newton's gravitational 
constant. The quantity G e g encodes the modification of gravity due to the presence of the 
scalaron field. Thus, the Poisson equation in Fourier space is transformed by replacing 
Newton's gravitational by the effective one in Eq. (4.29). The transition from the GR 
regime to the scalar-tensor one occurs when m ~ (aH/k) 2 . For the wave numbers relevant 
to the observable linear region of the matter power spectrum we require m(z = 0) > 3x 10 -6 
for the transition to have occurred by today. The Starobinsky model in particular allows 
for a rapid growth of m from R 3> Hq (m < 10~ 15 ) to R ~ Hq (m = O(0. 1)). 

Another way to understand the evolution of 7 in the Starobinsky model is analysing 
the evolution equation [48]. At high redshifts, we can approximate this evolution equation 



to [49] 




|d-7) 



( 



Ccff 

G 




|(2 7 -1)(«^ E + 1) 
(4.30) 



At early stages, the first term of the right hand side of Eq. (4.30) dominates. This is why 
7 decreases as G e g/G increases, even becoming negative. As we approach the present era, 
the second term starts to dominate making 7 increase again. 

In the case of the DGP model, the effect of the extra dimension affects both the friction 
term in the evolution equation for 5 m and the source term as, similarly to f(R) gravity, we 
can define an effective gravitational constant as 



where here (3 is the same as the parameter described in Eq. (4.2). As we can see in Fig. (lb), 
the parameter Qm is smaller for the DGP model than for the standard ACDM one. To 
compensate for this lack of matter density, the growth index is higher than for the rest of 
the models, as it is clear from Fig. (Id). 

In the case of the LTB model, we see in Fig. (lc) that the equation of state parameter w 
differs significantly from all the other models. The reason for this can be seen by looking at 
the formula used to calculate w(a), Eq. (4.27). The expansion rate is similar to the ACDM 
case since the parameters chosen for the void model considered have been chosen to fit 
the supernovae data. As we can see, the growth index for this model is not significantly 
different from the flat universe scenario and also tends to an asymptotic value given by 
7 = 0.573, very close to the predicted one in the ACDM case. 

As regards the wCDM model studied, we have used the values of wo = —0.9 and 
w a = 0.2 to calculate the growth index using Eq. (3.13). We have chosen these values to 
see how a considerably different model from the standard ACDM one can have a growth 
index that is barely distinguishable from it, as it can be seen in Fig. (Id). 

4.4.1 Ansatz for the growth index 

In this subsection we compare the different models' growth indices' redshift evolution with 
the Ansatz [71, 73] 



We have plotted in Fig. 2 the known redshift dependence of the growth index for the 
different models together with the corresponding values derived with the above Ansatz. As 
can be seen very clearly, for most models this is an excellent approximation. However, for 
f(R) models it fails to work at relatively nearby redshifts. For this reason we have adapted 
our Fisher matrix analysis for this specific case, see Section 6. 

5 Fisher matrix analysis prospects for a survey like Euclid 

The Fisher information matrix can be used to quantify the amount of information that an 
observable random variable X carries about a set of unknown parameters 9. 




(4.31) 



7 (a) = 7o +7a(l - a) . 



(4.32) 
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Figure 2: The redshift evolution of the growth index compared to the linear approximation 
in Eq. (4.32) (continuous lines), for all four classes of models. Note that while for most 
models the Ansatz gives an excellent description, for f(R) models it deviates significantly 
even for moderate redshifts. 

More specifically, let us assume that we have a data set arranged as a n-dimensional 
vector X = {x±, X2, ...x n } which can be thought of as a random variable with its probability 
distribution function L (X; 0) which depends on the model parameters 9 in a known way. 
The Fisher matrix is defined as the second derivatives of the logarithm of the probability 
distribution function with respect the parameters: 



where C = — InL. For a sufficiently large sample (in cosmology, a large survey), the central 
limit theorem ensures that the probability distribution function of the random variable X 
becomes Gaussian centered about their true value. In this case, in Eq. (5.1), the logarithm 
of the likelihood function L (X; 6) becomes a quadratic function of the variable X about 
its maximum. So, only under the assumption of a large sample, the inverse of the Fisher 
matrix becomes the covariance matrix and it can be used to estimate the errors on the 
parameters and furthermore to look at degeneracies between them. Then, with the help of 
the Cramer-Rao inequality, one can estimate the minimum variance of the parameter from 
the data set: the root square of the diagonal elements of the inverse of the Fisher matrix 
give the best errors on the parameters, ie. A9i = ••• 

This method is extremely useful, for instance, when we want to optimize the design of 
a survey, in the sense that it is possible to model such a survey in order to maximize its 
information content on the particular parameters we want to measure. 

5.1 Complete galaxy power spectrum 

There are two ways we can get information about dark energy parameters by looking 
at the galaxy power spectrum: first, by marginalizing over the main observables: the 
Hubble parameter, the angular diameter distance, the growth factor and the redshift space 
distortion parameter and projecting on to the dark energy parameter space; secondly, by 
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taking the full information of the observed galaxy power spectrum and deduce from the 
latter the parameters of interest. 

Following [50-52] we write the observed galaxy power spectrum as: 



D\(z)H r (z) 



P ohs (z;k,fi) = J r \' [ G 2 {z)b{z) 2 {l + PflY Por(k) + P shot (z) (5.2) 



where the subscript r refers to the values assumed for the reference cosmological model, i.e. 
the model at which we evaluate the Fisher matrix, b(z) is the matter bias factor, G{z) is the 
growth factor of matter perturbations. Also here -P s hot is the shot noise due to discreteness 
in the survey, [i is the direction cosine within the survey, Pq t is the present matter power 
spectrum for the fiducial (reference) cosmology. For the linear matter power spectrum we 
adopt the CAMB output [53]. 

The wavenumber k is also to be transformed between the fiducial cosmology and the 
general one (see [50-52] for more details). 

The distortion induced by redshift can be expressed in terms of the f3(z) factor which 
is related to the bias factor via: 

The (l + // 2 /?) 2 factor accounts only for linear distortion in the redshift space and it should 
be considered as a first approximation. More work beyond Kaiser's small-angle and the 
Gaussian approximation is needed, as discussed in [54-57]. 

The galaxy over-densities are assumed to trace the underlying matter distribution 
through a fraction called bias factor, b(z,k). This quantity could be arbitrary, it could 
even depend on both time and scale, see [58-60]. Usually it is assumed that the bias on 
large scales is independent on scale [61], hence in the matter power spectrum, this term 
appears as a multiplicative factor which modulates the overall amplitude of the galaxy 
power spectrum. We will assume here a Gaussian linear bias, with redshift dependence 
b(z) = vl + z, because it provides a good fit to the H a line galaxies in the near-infrared, 
which are the target of a Euclid- like survey, see [62] . 

The factor G{z) is related to the growth rate by f(a) = — 1 and it is usually 

assumed to be independent on scale: the late-time change in the expansion rate affects all 
scales equally. However, the growth factor may depend on the scale k, for instance allowing 
perturbations also in the dark energy sector. In the last case, the existence of the dark 
energy sound horizon will introduce a k dependence on the growth factor, see for instance 
[19, 20, 63]. However, here we assume a scale independent G(z) as we have neglected the 
dark energy perturbations. 

The total galaxy power spectrum including the errors on the redshift can be written 
as [52] 

P(z,k) = P ohs (z,k)e k2 ^ (5.4) 

where a r = 5z/ H{z) is the absolute error on the measurement of the distance and 8z is 
the absolute error on redshift. 

The Fisher matrix provides a useful method for evaluating the marginalised errors 
on cosmological parameters. Assuming the likelihood function to be Gaussian, the Fisher 
matrix is [64, 65] 

f k ™ dlogP(k) dlog P(k) k 2 
" l in ~Mi ^'tof** (5 - 5) 
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Parameters 


P(k) 


mP(k) WL 




1 


total matter density 


O h 2 


O h 2 


J 'mo ' <■ 


2 


total baryon density 






n b0 h 2 


3 


optical thickness 


r 


T 


T 


4 


spectral index 


n s 


n s 


n s 


5 


matter density today 




^mo 


^mo 


6 


equation of state parameter 


W 




Wo 


7 


equation of state parameter 


W\ 




W\ 


8 


rms fluctuations 








For each redshift bin 


9 


growth index 


7(2) or {70, 7a} 




7(2) or 70 


10 


Hubble parameter 




log 5(2) 




11 


Angular diameter distance 




log D A (z) 




12 


Growth factor 




logG(z) 




13 


z-distortion 




log/3(z) 




14 


shot noise 


Ps 


P s 





Table 1: Cosmological parameters for the complete and marginalized (mP(k)) galaxy 
power spectrum and weak lensing. 



where the #'s are the parameters shown in Tab. 1; the derivatives are evaluated at the 
parameter values of the fiducial model and V e g is the effective volume of the survey, given 
by: 

n(r)P(k,fx) \ 2 _ ( nP(k,n) 
n(r)P(k,n) + l) ~ \nP(k,fi) + 1, 



Here fx = k ■ r"/k, r is the unit vector along the line of sight and k the wave vector, and 
the last equality holds for an average comoving number density n. The highest frequency, 
kma.x(z), is evaluated at z of the corresponding bin and it is chosen so as to avoid non- 
linearity problems both in the spectrum and in the bias; we choose values from 0.11 /i/Mpc 
for small z bins to 0.25 fo/Mpc for the highest redshift bins. 

Any submatrix of F^ 1 gives the correlation matrix for the parameters corresponding to 
rows and columns on that submatrix. The eigenvectors and eigenvalues of this correlation 
matrix give the orientation and the size of the semiaxes of the ellipsoid confidence region. 
This automatically marginalizes over the remaining parameters. The square root of the 
diagonal elements will give the errors of the corresponding parameter. 

5.2 Marginalized galaxy power spectrum 

This approach consists on evaluating the Fisher matrix Fy first for the main observables 
shown in Tab. 1 and then project into the dark energy parameters, Wq, w a and 7. Since 
we want to propagate the errors to the cosmological parameters above, we need to change 
parameter space. This will be done taking the inverse of the Fisher matrix F^ 1 and 
then extracting a submatrix, called F~^ containing only the rows and columns with the 
parameters that depend on wq, w a and 7, namely H(z), Da(z), G(z) and (3(z). Then 
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we contract the inverse of the submatrix with the new set of parameters; the new Fisher 
matrix will be given by 

_ dp m „ dp n 
blj ~ % tmn dq t ■ [bJ) 

The square root of the diagonal elements of the matrix gives the errors on the parameters 
wo, w a and j(z). The growth index parameter will be discussed in more detail in the next 
section. 

5.3 Weak lensing 

Following [66, 67], the lensing potential is $l = ^+3?, which describes the deviation of light 
rays. As previously mentioned, we limit ourselves to scalar perturbations at linear order 
in the Newtonian gauge Eq. (3.1); as we are considering dark energy only at background 
level, the gravitational potential is simply given by Eq. (3.6). As $ = ^ it follows that the 
weak lensing potential is: 

k 2 $> L = 2^f^A m . (5.8) 

The convergence weak lensing power spectrum (which in the linear regime is equal to the 
ellipticity power spectrum) is a linear function of the matter power spectrum convoluted 
with the lensing properties of space. For a ACDM cosmology it can be written as [68] 



rco j r / tt p \ 



(5.9) 



where the W%s are the window functions, P n i [Pi (k, z)] is the non linear power spectrum at 
redshift z obtained correcting the linear matter power spectrum P; (k, z), see [18] for more 
details. 

The Fisher matrix for weak lensing is given by: 

Fop = / sky £ ( 2e+ ^ M d ( Py ) Cj^d {P km ) tP C~] (5.10) 



where the partial derivatives represent d/d8 a , the corresponding cosmological parameters 
9 a are shown in Tab. 1 and 

C jk = P jk + 6 jk ^^ (5.11) 



1/2 

where 7^ is the rms intrinsic shear (here we assume (7 int )=0.22 [69]) and nj is the number 
of galaxies per steradians belonging to the i-th bin. 



6 Results 



We are now in the position to derive the sensitivity of a typical next-generation survey to 
the dark energy parameters. 

In particular, we consider a survey patterned according to the specification of the Eu- 
clid survey, see [13]. In order to do this we make use of the Fisher matrix formalism for the 
three different observables introduced in the previous section: complete and marginalized 
galaxy power spectrum and Weak Lensing. 
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z 


ni(z) x 10 3 


n 2 (z) x 10 3 




0.5-0.7 


4.69 


3.56 


0.7-0.9 


3.33 


2.42 


0.9 - 1.1 


2.57 


1.81 


1.1 - 1.3 


2.1 


1.44 


1.3 - 1.5 


1.52 


0.99 


1.5-1.7 


0.92 


0.55 


1.7-1.9 


0.54 


0.29 


1.9 - 2.1 


0.31 


0.15 



Table 2: Expected galaxy number density in each redshift bin for the Euclid survey in 
units of (/i/Mpc) 3 . 



Galaxy redshift survey specifications: for the galaxy power spectrum cases we 
consider a spectroscopic survey from z = 0.5 — 2.1 divided in equally spaced bins of width 
Az = 0.2 and with a covering area of 20000 deg 2 ; the galaxy number densities in each bin 
are shown in Tab. 2 with an efficiency of 50%. For all cases we assume that the error on 
the measured redshift is spectroscopic, 5z = 0.001 (1 + z). 

While modeling the redshift survey, we choose two different cases according to the 
galaxy number density [70, 71]: 

• optimistic case; this corresponds to the middle column in Tab. 2. 

• realistic case; this corresponds to the last column in Tab. 2. 



Weak lensing survey specifications: for the WL survey we consider a photometric 
survey characterized by the sky fraction / s k y = 1/2, that is, a covering area of 20000deg 2 ; 



with zq = Zmean/1-412 and mean 
axies per arcmin 2 . Moreover, we 



an overall radial distribution n(z) = z 2 exp — (z/zq) 1 ' 5 

redshift z me an = 0.9, the number density is d = 35 ga... 
consider the range 10 < I < 10000 and we extend our survey up to z max = 3 divided in 
5 bins each containing the same number of galaxies. For the non linear correction we use 
the halo fit model by Smith et al. [72]. We assume the error on the measured redshift is 
photometric, 5z = 0.05 (1 + z). 

Fiducial model: our fiducial model corresponds to the ACDM WMAP-7yr [6]: 
n m>0 h 2 = 0.134, n b h 2 = 0.022, n s = 0.96, r = 0.085, h = 0.7, ft m , = 0.275 and U K = 0. 
For the dark energy parameters we choose Wq = — 1 and w a = 0. 



We can now derive the sensitivity of the parameters introduced above. More specif- 
ically, we consider separately three different cases: in Case 1 the growth index has been 
chosen to be a free parameter but independent in each redshift bin; in Case 2 we consider 
the growth index as a free parameter but equal for all the redshift bins; in Case 3 the 
growth index depends directly on the equation of state parameters (wq and w a ), using 
Eq. (3.18) as the analytic solution. 
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Figure 3: Marginalized errors for the growth index 7 at different redshift bins for the 
complete (left panel) and marginalized (right panel) galaxy power spectrum. The blue 
dashed and the red errors bars refer to the optimistic and realistic case, respectively. The 
black long dashed line is the growth index give by the Eq. (3.18). 




Figure 4: Marginalized errors for the growth index 7 at different redshift bins for the 
WL convergence spectrum. 



• Case 1: 

We consider the growth index 7 as a free parameter and independent for each redshift 
bin in order to map its variation over time; its value is chosen to be 0.545. The errors 
are shown in Fig. 3 for both galaxy power spectrum cases; in Tab. 3 are reported 
their 1 — a errors. Overall the errors on the growth index are of about 0.02 and 0.03 
for the complete and marginalized galaxy power spectrum, respectively. In Fig. 4 the 
errors for the WL case are shown. It is worth noticing that after z ~ 1 the errors on 
7 are basically unchanged. Furthermore, we also plot the analytic expression for the 
growth index Eq. (3.18) as a comparison. However, the difference between the full 
analytic expression for the growth index and its asymptotic value 7 = 0.545 is too 
small to be detected even for a half sky survey like Euclid for both the spectroscopic 
and the photometric cases. 

• Case 2: 
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P(k) 


z 


0.6 


0.8 


1.0 


1.2 


1.4 


1.6 


1.8 


2.0 


n real 
7 


0.0305 


0.0231 


0.0194 


0.0159 


0.0153 


0.0162 


0.0197 


0.0281 


opt 


0.0296 


0.0219 


0.0183 


0.0147 


0.0139 


0.0141 


0.0153 


0.0185 


mP(k) 


z 


0.6 


0.8 


1.0 


1.2 


1.4 


1.6 


1.8 


2.0 


n real 

7 


0.0311 


0.0286 


0.0284 


0.0308 


0.0337 


0.0367 


0.0397 


0.0427 


opt 
(T-y 


0.0286 


0.0256 


0.0250 


0.0268 


0.0292 


0.0318 


0.0342 


0.0366 



Table 3: Here are listed the 1 — a errors for the growth index 7 at different redshifts for 
the P (k) methods. 



WL 


z 


0.5 


1.0 


2.0 


3.0 




0.097 


0.055 


0.051 


0.051 



Table 4: Here are listed the 1 — a errors for the growth index 7 at different redshifts for 
the WL method. 

For this case we chose the growth index to be [71, 73] 

70) = 7o +7o(l - a) = 7o + 7a— — (6.1) 

1 + z 

in order to compare the alternative theories introduced in the previous sections. Our 
reference cosmology is always ACDM with 70 = 0.55620 and 7 a = 0.0182537. We 
choose these values according to the Taylor's expansion of Eq. (3.18) for the ACDM 
cosmology. In Fig. 5 the 1 — a confidence regions for the parameters wq — w a (left 
panel) and 70 — 7 a (right panel) for the galaxy power spectrum cases are shown. 

For the WL case we found j a to be degenerate with 70 indicating probably that 
a weak lensing experiment is insensible to the variation of the growth index with 
redshift; this can be also seen in Fig. 4 where the errors of the growth index are 
basically the same for all the redshift bins above z — 1. In Tab. 5 we report the 
errors for the parameters introduced above for all the three surveys. 

We reach a sensitivity of about a WQ = 0.018, a Wo = 0.068 and a WQ = 0.122 for 
the galaxy power spectrum cases and WL survey, respectively. These errors are 
sufficiently small to rule out independently most of the models we considered in this 
paper, see Fig. 5 where we also plot the expected values of the equation of state 
parameters wq — w a and the growth index 70 — 7 a for the alternative cosmological 
models: quintessence model (yellow box), DGP (blue diamond), LTB (black triangle) 
and f{R) model (brown inverted triangle). Only a measurement of the equation of 
state parameter is able to rule out wCDM, DGP and LTB model; however, f(R) 
is almost indistinguishable from the reference cosmology and none of the surveys 
assumed here are able to rule it out. Fortunately it is also possible to measure the 
growth index 70; in this case we reach a sensitivity of about 0.02, 0.092 and 0.075 for 
the galaxy power spectrum cases and WL survey, respectively. On the other hand, 
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Figure 5: Confidence level for the wq, w a (left panel) and 70 , j a (right panel). The 
blue and light blue areas are for the galaxy power spectrum case, optimistic and realistic 
case respectively. The red dashed and solid lines are for the mP(k) case, optimistic and 
realistic case respectively. The light green shaded area is for the WL case. These two plots 
correspond to case 2 of our analysis, taking 7 = 70 + 7a(l — a). Note that the vertical 
axis has been broken on the left panel plot to include the point corresponding to the LTB 
model. The points corresponding to the different models are the wo,w a and 70, 7a pairs 
calculated for each model. 





P(k) 


mP(k) 


WL 




real. 


opt. 


real. 


opt. 






0.021 


0.018 


0.076 


0.068 


0.122 


°Wa 


0.051 


0.041 


0.375 


0.324 


0.524 


<T 7o 


0.022 


0.020 


0.102 


0.092 


0.075 


°7a 


0.120 


0.116 


0.339 


0.296 





Table 5: Here are listed the 1 — a errors for wq, w a , 70 and 7 a for the P (k), mP(k) and 
WL cases. 

measuring only 70 — 7 a , the model that could be clearly ruled out is DGP, while LTB 
could still be a viable model. However, since we are dealing with two different surveys 
we can consider measuring one set of parameters for each experiment; for example, 
we can measure wq — w a with the WL survey and 70 — 7 a with a P(k) experiment. In 
this case the models considered in this work would be excluded by at least at the 2a 
level. With respect to f(R) models, note that the Taylor expansion of 7(a) does not 
work appropriately in the whole range of redshifts explored by the survey. Therefore, 
we have left this class of models out of Fig. 5. However, inspecting Fig. 7 one can 
easily conclude that models f(R) would also be ruled out by the measurements of a 
Euclid-like survey at the many a level. 

• Case 3: 

In this case, the growth index 7 depends on the cosmological parameters wq — w a 
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Figure 6: Confidence level for wo and w a for all the three cases with the growth index 
given by Eq. (3.18). The red solid line and blue dashed line refer to the optimistic and 
realistic case, respectively. Note that the vertical axis has been broken on the left panel plot 
to include the point corresponding to the LTB model. Note that the right figure contains 
a scaled version of the left figure. For details see the text. The points corresponding to the 
different models are the Wq, w a and 70, 7a pairs calculated for each model. 



according to Eq. (3.18); the derivatives of the growth factor are given by: 



dlogG 



9r, _ , , d log Q(z) 

io g n m (z) +7- 



-'"to - "III V, / 1 ' o 

dqi dqi 



(6-2) 



which has an extra term dj/dqi, our new set of parameters becomes now qi = 
{wo,w a }- In Fig. 6 we plot the 1 — a confidence regions for wq — w a for all the 
three experiments; in the same figures are also plotted the values of the equation of 
state parameters WQ — w a for the alternative cosmological models: quintessence model 
(yellow box), DGP (blue diamond) and LTB model (black triangle) 2 . 

However, in order to be able to include in the same graphic the other cosmological 
models we had to multiply the confidence regions by a factor of 5 for the galaxy power 
spectrum cases and by a factor of 3 the WL case; in the right panel of Fig. 6 we zoomed 
over the confidence region to preserve the proportions of the errors obtained with the 
three different methods. Giving a direct dependence to the growth index 7 on the 
equation of state parameters the errors on wo — w a are reduced by a factor 4 for the 
P{k) case, a factor 10 for the mP(k) and of about a factor 2 for WL. In this case we 
reach an extreme sensitivity of about 0.5 — 1.0% for wo and of about 2 — 3% for w a , 
however still not sufficiently good to rule out completely the f(R) model studied. 



2 We have not plotted the fs (R) model here because the Ansatz for 7(a) is not valid in the whole range 
of redshifts. 
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P(k) 


mP(k) 


WL 




real. 


opt. 


real. 


opt. 






0.0052 


0.0047 


0.0063 


0.0056 


0.072 




0.0155 


0.0135 


0.0281 


0.025 


0.106 



Table 6: Here are listed the 1 — a errors for wq and w a for the P(k), mP(k) and WL 
cases. 

7 Conclusions 

Over the last decade there has been a plethora of proposals to account for the observed 
acceleration of the universe. All these proposals fall into four main categories: i) the 
inclusion of some extra field (scalar, vector or tensor), coupled or not to the rest of matter, 
like in quintessence, chameleon, vector dark energy or massive gravity; ii) the extension 
of GR by inclusion of higher order terms in the Einsten-Hilbert action, like f(R) theories, 
Gauss-Bonnet terms, etc.; iii) the modification of gravity on large scales by introduction of 
extra dimensions, like in Dvali-Gabadadze-Porrati model, Kaluza-Klein gravity, etc.; iv) the 
reinterpretation in terms of a nontrivial spatial geometry, like in large- void inhomogeneous 
LTB models. 

All of these proposals have very specific predictions for the background evolution of 
cosmological space-times, and most of them can be well fitted to the present observations, 
with just a few phenomenological parameters: the equation of state, the speed of sound, 
the coupling between DM and DE, bulk viscosity, etc. However, in order to discriminate 
between the different alternatives it has been realized that one has to go beyond the back- 
ground evolution and start to consider also the theory of linear cosmological perturbations 
and parametrize their evolution in terms of the growth function and growth index, as well 
as the shift parameter. 

In this paper we provide exact solutions to the cosmological matter perturbation equa- 
tion in a homogeneous FLRW universe with a vacuum energy that can be parametrized with 
a variable equation of state parameter w(a) = wo + w a (l — a). We compute the growth in- 
dex 7 = log /(a) / log U m (a), and its redshift dependence, using the exact solutions in terms 
of Legendre polynomials and show that it can be parametrized as 7(a) = 70 + 7 a (l — a) 
for most cases, see Fig. 2. We then compare four different types of dark energy models: 
iuACDM, DGP, f(R) and a LTB-large-void model, which have very different behaviors at 
z>l, see Fig. 7. This allows us to study the possibility to differentiate between various 
alternatives using full sky deep surveys like Euclid, which will measure both photometric 
and spectroscopic redshifts for several million galaxies up to redshift z = 2. We do a Fisher 
matrix analysis for the prospects of differentiating among the different DE models in terms 
of the growth index, taken as a given function of redshift or with a principal component 
analysis, with a value for each redshift bin, see also [74, 75] for a similar analysis for DGP 
and f(R) models. We use as observables the complete and marginalized power spectrum of 
galaxies P(k) and the Weak Lensing (WL) power spectrum. We find that using P(k) one 
can reach (2%, 5%) errors in (wo,w a ) and (4%, 12%) errors in (70,7a), while using WL we 
get errors at least twice as large. These estimates allow us to differentiate easily between 
DGP, f(R) Starobinsky model and ACDM, see Fig. 7, while it would be more difficult to 
distinguish the latter from a variable wCDM or LTB models using only the growth index. 
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Figure 7: The redshift evolution of the growth index compared with the expected errors 
from a Euclid-like survey. Note that DGP and the studied f(R) models can in principle 
be ruled out at many sigma. 
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Appendices 

A Small scale approximation 

In this appendix we present the small scale solution found to the perturbation equation 
for the matter density contrast. The master equation for the perturbations in the small 
scale regime (i.e.: large wave number k) (3.7) can be solved by changing the independent 
variable from time to the scale factor, 



a 2 5"(a) + (3 - e(o)) a 5' (a) - ~Q m {a) 5(a) = . 



where e is given by 



e(a) 



aH' _ 3 
'IT ~ 2 



1 + w(a)[ 1 - TO (o) 



(A.l) 
(A.2) 
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Furthermore, by noting that the density contrast for pure matter grows like 5(a) oc a, one 
can write the general function 

5(a) = a • G(a) , 



with which the master equation becomes 

a 2 G"(a) + (5 - e(o)) a G'(a) + (3 - e(o) - ^ m (a)) G(a) = . 
Now we make a change of variables, 

^de,0 



-3w 



m,0 



and thus the master equation becomes 



u(l-u)G"(u) + 



5 

Qw 



da 



-3w u- 



-) 



u 



du 



G'(u) 



n m (o) 



w 



Qw 2 



G(u) = . 



(A.3) 



(A.4) 



(A.5) 



(A.6) 



which has the form of a Hypergeometric equation, see Ref. [23], with constant coefficients 
a = (w — l)/2w,f3 = — l/3u;,7 = 1 — 5/6w, and thus the exact solutions are written in 
terms of two independent constants, C\ and C2, 



5(a) = Ci a ■ 2F1 



w 



1 -1 



2w ' 3u; 
1 1 



,1 



5 _ 



^de,0 



-3w 



m.O 



(A.7) 



zu; 2 3w; 



-3ui 



The first term corresponds to the growing mode solution and the second one to the decaying 
mode. When describing late time solutions we will always take the growing mode solution; 
furthermore, the integration constant C\ is not a problem, as in this work we are interested 
in the evolution of the growth rate f(a) which is the ratio of the matter density and its 
derivative. 

There are a number of way that one can play with the solution above in order to simplify 
the Hypergeometric function and even to drop the Hypergeometric functions. Here we use 
one relation which seems to be the easiest among all the others; we can notice that the 
third coefficient 7 of the Hypergeometric function can be written as: 



1 



7 



+ a + /3 



(Ai 



being a and /3 the first and second coefficient of the Hypergeometric function. In this case 
the solution for the matter density contrast can be written in terms of Legendre polyno- 
mials. However, being our goal to find an exact solution for the growth rate and growth 
index, we find easier to first evaluate these terms using the Hypergeometric functions and 
then to simplify the result. The growth rate is defined as 



/(a) 



a5'(a) 
~J(aJ 



1 + a- 



af3n> m (a) 2*1 


a + l,/3 + 1,7 + l,l-n- 1 (o) 


7 ^m(a) 2Fl 


a,/3,7, 1 - Vim (a) 





(A.9) 
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using Eqs. (15.4.12) and (15.4.21) of [23] we have: 
f{a) = 1 + 6wa(3 y/l - fi m ( 



l/v/«m(a) 



pi 



(A.10) 



being P™ (x) the Legendre polynomial. The last equation can be further simplified making 
using of the recurrence relations of the Legendre polynomials, see Eqs. (8.5.3) and (8.5.5) 
of Ref. [23] , then we find 



/(a) = nU\a) 



p5/6«« 
l/6w 


O m 1/2 (a) 


p5/6to 

-1/&W 


^m 1/2 (a) 



(A.ll) 



We can then express the growth index 7(a) in terms of the Legendre polynomials, and we 
find: 



7(a) 



1 



+ 



1 



In fi m (a) 



In 



p5/6to 
l/6w 


^m 1/2 (a) 


p5/6«i 

-1/&U) 


^m 1/2 (a) 



(A.12) 



Note that the growth index parameter is close to 0.5 ± 0.1, depending on the range of 
values of the equation of state parameter and the matter content of the universe. 

A.l Varying equation of state parameter 

We can extend our discussion also to a varying dark energy equation of state parameter 
w(a). It is fair to be said that in the last case there is no exact analytic solution for the 
matter density contrast as it was for the case in which w is constant; the main problem 
here is that there is no a direct transformation between the scale factor a and the new 
variable u. 

Here we assume the equation of state parameter to be: 



w(a) = w + w a (1 - a) , 
for which the matter density parameter can be integrated 



fi m (a) 



de,0 n -3(wo+Wa) p 3w a (a-l) 



-1 



(A.13) 



(A.14) 



The master equation for a varying equation of state parameter still looks the same as 
Eq. (A. 4), except that now w is a function of the scale factor. In order to integrate out we 
use again the variable 



1 - O 



u 



n 



nx,0 a -3w(a) 



m,0 



with w(a) 



1 

In a 



a w (a r ) 



da . 



(A.15) 



However, we can make the approximation that the w(a) is slowly varying with time and 
integrate out the master equation. The approximation is not so rude. To see this, we need 
to have another look at the master equation for matter perturbations. 
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In this case the dark matter density contrast reads 



5(a) 




(A.16) 



with parameters 



w(a) — 1 w a a 



(A.17) 



a 



2w(a) 6w 2 (a) 
-1 w a a 



3w(a) 6w 2 (a) 



(A.18) 



which gives a density growth function 



/(a) = nU\a) 



■,5/6w(a)+w a a/6w 2 (a) q— 1/2 
1/Qw{a) [ ilm 



(A.19) 



) 5/Gw(a)+w a a/6w 2 (a) q— 1/2 
-l/6u>(a) |_ SZm 



Comparison of the numerical solution for the density growth function with the approximate 
expression (A.19) shows accordance within less than 0.1% for a very wide range of values 
of -2 < w a < 2. 
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